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Abstract. We calculate the tidal interaction of a uniformly rotating 1 Mq star with an orbiting companion at 
various phases of core hydrogen burning from the ZAMS to core hydrogen exhaustion. By using the traditional 
approximation we reduce the solution of the non-adiabatic oscillation equations for the tidal forcing of a rotating 
star to a one dimensional problem by solving a separate eigenvalue problem for the angular dependence of the 
tidal perturbations. The radial oscillation equations are then solved by using finite differencing on a fine grid 
with an implicit matrix inversion method like for stellar evolution calculations. We are able to identify resonances 
with gravity (g-)modes and quasi-toroidal (q-)modes with up to ~ 1000 radial nodes in the more evolved stellar 
models. The resulting tidal torque is calculated down to low forcing frequencies close to corotation. For low 
tidal frequencies we find significant response due to inertial (i-)modes in the convective envelope. The inertial 
modes are damped by turbulent dissipation in the envelope and generate a relatively high torque-level in the low 
frequency region where the (retrograde) high radial order g-mode resonances become tidally inefficient due to their 
rotational confinement to the stellar equator and strong damping by radiative losses. For still lower retrograde 
forcing frequencies we find a large number of closely spaced weakly damped quasi-toroidal q-mode resonances. 
Our results indicate that effects related to stellar rotation can considerably enhance the speed of tidal evolution 
in low-mass binary systems. 

Key words. Hydrodynamics- Stars: binaries- Stars: rotation- Stars: oscillation 



r-j . 1. Introduction 



Spectroscopic observations show there is a fairly well de- 
termined orbital period below which main sequence bina- 
ries with solar-like stars have essentially circularized or- 
bits. This circularization period appears to be P c i rc ~ 11— 
12 days (Duqucnnoy and Mayor 1991; Latham ct al. 1992; 
Mathieu et al. 1992| ). In wider binary systems the tidal 



eluded that it predicts a circularization period of only 
P c i rc ~ 2.2 days if one allows for the reduced efficiency 
of the viscous damping when the convective timescale of 
turbulent eddies is longer than the typical tidal forcing 
period. They find P c irc — 6 days if one ignores this effect. 
Because the tidal torque scales as P~£, where P or b is the 
orbital period, this forms a significant discrepancy with 
the observationally inferred P c ; rc ~ 11 days. 



evolution is believed to be too slow to achieve circulariza- 
tion during the main sequence phase of solar- type stars. 
Since Zahn's (1977) review of tidal effects there has re- 
cently been a renewed interest in tides in solar-type bi- 
nary system, because the circularization times predicted 
by Zahn's analysis are too long to explain the observa- 
tions. Zahn considered the interaction of turbulent ed- 
dies with the hydrostatic equilibrium tide, i.e. the low 
frequency limit of the tidal response, as the dominant 
tidal mechanism in low mass stars with convective en- 
velopes. Although this mechanism seems consist ent with 
observed binary sy stems containing giant stars ( Vcrbunt 
and Phinney 199£ ) , it appears to fall short by orders of 
magnitude when co mpared with observed so l ar-type bi- 
nary systems (e.g. Claret and Cunha 1997 ). |Goodman 



and Oh (1997 ) also applied Zahn's formalism and con- 



Tcrquem et al. (1998, from now on T9S) calculated 
the dynamical tide raised by a close companion on a non- 
rotating solar-type star. For dynamical tides the interac- 
tion between the tidal flow and the stellar normal modes 
of oscillation is taken into account flCowling 1941 ). The 
stratification in low mass main sequence stars is comple- 
mentary to that of massive stars and consists of a ra- 
diative core surrounded by a convective envelope. Near 
the boundary of the convective envelope the tides can ex- 
cite g-modes in the radiative core, where radiative damp- 
ing is relatively weak. Additional damping occurs in the 
convective envelope by the interaction of the tidal flow 
with turbulent eddies. T98] allowed for damping due to 
the interaction with these turbulent convective eddies by 
first order perturbation analysis, while radiative damping 
was treated by a WKB treatment of non-adiabatic terms. 
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They noted that the equilibrium tide applies to forcing fre- 
quencies much smaller than the Brunt- Vaisala frequency 
|7V 2 | 5 (see below) and that for solar-type stars with rel- 
atively small values of | AT 2 | 5 in the lower regions of the 
convective envelope the convective timescale is too long 
for essentially instantaneous adjustment to the tidal flow. 
Under these circumstances the equilibrium tide is a rather 
poor approximation for the tidal interaction and predicts 
a too strong effect. T98j found that for a fixed spectrum 
of normal modes in a non-rotating star the tidal evolu- 
tion of solar-type binary systems is controlled essentially 
by the non-resonant dynamical tide. To explain the ob- 
served P c irc would require a viscosity that is ~ 50 times 
larger than predicted by sim ple mixing length estimates. 
Goodman and Dickson (1998 ), in their WKB treatment of 



dynamical tides, considered non-linear damping of excited 
g-modes near the stellar centre but they estimate this can- 
not explain the discrepancy with the observations. 



In T98 the stellar normal mode spectrum was consid- 
ered fixed, i.e. for simplicity effects of stellar rotation and 
stellar evolution were ignored. However, in orbits with a 
significant eccentricity the tidal forcing is characterised 
not by one but by a spectrum of forcing frequencies which 
can potentially excite one or more of the normal modes 
in the binary stars. It is often argued that the chance 
of hitting a (narrow) resonance is small and that -when 
it happens- the enhanced tidal interaction rapidly moves 
the system through the resonance without much orbital 
evolution. This ignores the possibility of resonance lock- 
ing. When stellar rotation is taken into account, the tides 
can excite both prograde (propagating in the direction of 
rotation) and retrograde (counter to rotation) oscillation 
modes, depending on the stellar rotation rate and orbital 
parameters. Excitation of prograde modes gives rise to 
spin- up (positive torques), while retrograde tidal forcing 
yields negative t orques and spin-down. Ca lc ulation s for 
massiv e binarie s ( Witte and Savonije 1999^ Jb|, |200l| , from 
now on |WS99a| ; |WS99b| ; |WS0l|) indicate that often a weak 
high order prograde orbital harmonic is driven into reso- 
nance with a stellar gravity (g-)mode by the simultaneous 
action of a strong low frequency retrograde orbital har- 
monic which spins the star down, increasing the relative 
forcing frequency of the high orbital harmonic towards res- 
onance. A similar driving effect can be caused by stellar 
evolution induced spin-down. The near-resonant weak pro- 
grade orbital harmonic tend to spin the star up, counter- 
acting the retrograde driving ever more strongly when the 
peak of the resonance is approached. The two counteract- 
ing (retrograde and prograde) orbital harmonics then tend 
to lock each other at resonance with prolonged, consider- 
ably enhanced, tidal interaction. The effectiveness of the 
locking process is sensitive to the magnitude of the orbital 
eccentricity and to the low-frequency retrograde tidal forc- 
ing. Unfortunately, the low-frequency forcing excites high 
order (strongly damped) retrograde g-modes with short 
wavelengths which requires a very large number of mesh- 
points for finite difference calculations. The 2D oscillation 



code (Savonije and Papaloizou 1997) which accounts for 
the full Coriolis force is not suited for the low frequency 
calculations with low-mass stars. 

In this paper we calculate the fully non-adiabatic tidal 
interaction with a uniformly rotating solar-type star by 
applying the traditional approximation for which the z?- 
component of the stellar angular velocity is neglected. This 
is a reasonable approximation for low oscillation frequen- 
cies which simplifies the tidal problem significantly in that 
it allows a full separation of the tidal perturbations in the 
three spherical coordinates, like for a non-rotating star. 
Here we present tidal torque calculations (down to low 
forcing frequencies) for a uniformly rotating 1 M star at 
various stages of core hydrogen burning. 

2. Basic tidal oscillation equations 

We consider a uniformly rotating main sequence star with 
mass M s = 1 Mq and radius R s . We assume the star's an- 
gular velocity of rotation H s to be much smaller than its 
break-up speed, i.e. {n s /fl c ) 2 < h with fl 2 = GM S /R%, 
so that effects of centrifugal distortion (oc il 2 ) may be 
neglected to a first approximation. We wish to study the 
response of this uniformly rotating star to a perturbing 
time-dependent tidal force. The Coriolis acceleration is 
proportional to O s and we consider its effect on the tidally 
induced motions in the star. We use nonrotating spheri- 
cal coordinates (r, (/?), with the origin at the primary 
star's centre, whereby $ = corresponds to its rotation 
axis which we assume to be parallel to the orbital angular 
momentum vector. 

Let us denote the displacement vector in the star 
by £ and perturbed Eulerian quantities like pressure P', 
density p' , temperature T" and energy flux F' with a 
prime. The linearised hydrodynamic equations governing 
the non-adiabatic response of the uniformly rotating star 
to the perturbing potential <!>t may then be written as 



1 , p' 
p p 1 



d_ 

dt 



El 
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dT 
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■ £ ■ vs] 
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(2) 

(3) 
(4) 



where n denotes the opacity of stellar material and S its 
specific entropy. These perturbation equations represent, 
respectively, conservation of momentum, conservation of 
mass and conservation of energy, while the last equation 
describes the radiative diffusion of the perturbe d energy 
flux. For simplicity we adopt the Cowling (1941 ) approxi- 
mation, i.e. we neglect perturbations to the gravitational 
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potential caused by the primary's distortion. We also ne- 
glect perturbations of the nuclear energy sources and of 
convection. 

Let us here consider the simplest case of a close binary 
system in which the orbit is circular with angular velocity 
oj and orbital separation D. The dominant tidal term of 
the perturbing potential is then given by the real part of 

$ T (r,^^t) = /r 2 P 2 2 ( M )e i ( CTt - 2 ^ 



where a — 2lo is the forcing frequency in the inertial frame, 
[i = cos i?, P|(m) i s the associated Legendre polynomial 
and / = — -Tjpr , with M p the companion's mass. After 
separating the ^-dependence, Eqs. (|l|-^) form a 2D prob- 
lem in the (r, i?) meridional plane of the perturbed star. 
Let us now apply the traditional approximation in order to 
attain a considerable simplification of the tidal problem. 

3. The traditional approximation 

For Sl s ^ the solutions of Eqs. arc no longer fully 

separable into r-, and ip- factors due to the Coriolis 
force. However, in the traditional approximation separa- 
bility is retained b y neglecting the component of the an- 
gular velocity (e.g. Unno ct al. 1989 ). This is done because 
the radial motions are expected to be small in the stably 
stratified layers of the star, especially for low-frequency 
high-order g-modes. For this reason, it is thought to be a 
reasonable approximation for low oscillation frequencies. 
Furthermore, a spectrum of rotationally governed inertial 
modes as well as (q uasi-)toroidal modes is also present in 
this approximation ( [Bavonije et al. 1995 ). 

When only the radial component of J7 S is retained, the 
$- and (^-components of the equation of motion can be 
written: 

1 dP' 



a — 2iSl s <TCOSi9^ 



a + 2iQ s <7cos$£ij 



rp dfi 



i 

~r~~d¥ 



rp sin I? 



• sin 



<j> n 



(6) 



(7) 



We have expressed the velocity perturbations in terms of 
the displacement vector by the relation v' = icr£, with 
a = a — mil s the oscillation frequency in the corotating 
frame, where m = 2 in our study. We assume m to be 
always positive, whereby a retrograde oscillation (propa- 
gation direction counter to the sense of stellar rotation) 
corresponds with negative oscillation frequencies a. 

3.1. Separation of variables 

A separation of variables can be performed by writing, see 
Papaloizou and Savonije (1997): 



n—1 n—1 

oo 

P' = J2Xn(#)W n (r), (8) 



with £ r , T 1 and p' having expansions of the same form as 
that for P' . Here T n , Q n and X n are functions of $ chosen 
to obey the relations 



-a T n — 2if2 s <r cos -dQn 



dXn 

' d-d 



and 



( 5 ) -a 2 G n + 2ifl s o- cos dT n 



lm 
sin 



X.„ 



We obtain an equation for X n (d) by imposing the con- 
straint 



1 d(sint?^" n ) . Q n \ n 

■ — 1H! = rrA ? , 



where A„ is a constant. Then X n must satisfy the second- 
order equation obtained from 



1 dQn , mx cosd _ I rrv 



sin$ d$ sin 2 i9 " \sin 2 $ 



with 



Qn 



sini? 



(1 - x 2 cos 2 i?) V di? 



dX, n mx cos $ 



sin 



X, 



(9) 



(10) 



whereby X n (■&) is an eigenfunction with X n the associated 
eigenvalue. The solution depends on the rotation parame- 
ter x = 2fl s /a. The tidal problem we are interested in cor- 
responds with even solutions of X n {p) with p = cosi? and 
m = 2. For Sl s = the functions X n {"&) become the asso- 
ciated Legendre functions P 2 2 „ (cos $) with corresponding 
eigenvalues A„ = n (n+ 1). Normal modes of the rotating 
star correspond to normal modes of the non-rotating star 
with n (n + 1) replaced by any permissible A„. 

Different X n ($) are orthogonal on integration with re- 
spect to p = cos$ over the interval (—1, 1). If the perturb- 
ing potential is expanded in terms of the X n such that 



oo 

$ T M) = 53Vn(r)X n (<>), 

n=l 

we find from (^J) and (Q) that 

D n (r) = —W n (r) + -V»(r). 
rp r 



(11) 



(12) 



n=l 



This equation is exactly the same as in the non-rotating 
case. The same is true for the equation of continuity, ex- 
cept that n(n + 1) is replaced by A n . If we ignore the en- 
ergy equation (see below) the tidal response will ensist of 
a superposition of responses appropriate to non-rotating 
stars with n (n + 1) replaced by A„. Because the angular 
functions X n ($) are not spherical harmonics for non-zero 
rotation, all n = / are required in the superposition. 

3.2. The modified oscillation equations 

Once we have solved the above eigenvalue problem for A„ 
the linearised equations which describe the non-adiabatic 
forced tidal oscillations for the n-th tidal component in 
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expansion ( [TT| ) become a 1-dimensional (radial) problem 
by factoring out the common ^-dependence factor X ra ($) 
of all occurring perturbations. The oscillation equations 
can now be expressed as 



dP' 
dr 



dP 
dr 



1 d(pr 2 £, r 
dr 



pr* 



P_ 

P 



P' p' 

P P 



2pfr, 



1 d(r 2 F;) 
Fr 2 dr 



+ fr 2 



(13) 



(14) 



,(15) 



3.3. Turbulent viscosity 

The tidal oscillations are expected to be influenced by the 
turbulent convective motions in the stellar envelope. We 
describe these interactions in a heuristic way by means of 
a simple mixing length approximation for turbulent vis- 
cosity. To this end we add in the convective envelope a 
radial viscous force density 



, _ icr d 

Jr. vise — ~~ o" ~7T~ 

r or 



pvr 



Mr 

dr 



to the equation of motion (|13|), with v a turbulent viscosity 
coefficient. In order to enable comparison we adopt the 
same form for the turbulent viscosity coefficient as T98: 



PL . d (V\ IA . T 1 ,„ .p' K X dI, ,,„, 
f = A dr [ T ) + ^ ] T-^- p + ^^ (16) 



where T's represent Chandrasekhar's adiabatic coeffi- 
cients, X the hydrogen abundance, 



A 



K T 



dlnT 
dr 



9hiK 
<91nT 



A 



dlnp 1 dlnP 
dr Ti dr ' 



d In k 
dh\p 



KX 



d\nn 



-(r 3 - 1)F/ {aP) , with F 



AacT 6 dT 
3np dr 



the unperturbed radiative energy flux, and where the 
other constants have their usual meaning. The factor r\ is 
a characteristic radiative diffusion length. Note that the 
radiative diffusion introduces a factor i, so that the tidal 
perturbations are complex-valued which expresses the in- 
duced phase-lags with respect to the companion. 

Although the energy equation is not separable, we have 
nevertheless simply replaced the n(n + 1) factor in the en- 
ergy equation for a non-rotating star by A„ (eq. |l^) to 
approximately take the horizontal radiative diffusion into 
account. It appeared that the inclusion (exact for non- 
rotating stars) or non-inclusion of this term yields almost 
the same values for the tidal torque, even when |A| be- 
comes large at low retrograde forcing frequencies. 

Finally, the oscillation equations are complemented by 
the linearised equation of state 



~P 



d\nP\ T 
dWT ) ~f 



ainP 

dlnp 



P_ 

P 



din p, a J p, a 



where p a is the mean atomic weight of the stellar gas. In 
radiative layers we assume convection of oscillatory stellar 
material without any diffusion, so that 



Pa 
Pa 



d In fig 
dr 



L 2 

7~con 



l 



( Tcon 27r) 



(18) 



where L is the mixing length, taken to be twice the lo- 
cal pressure scale height, r con = |A/" 2 |~2 the characteris- 
tic convective timescale with M 2 = —gA the square of 
the Brunt- Vaisala frequency, whereby g denotes the local 
acceleration of gravity. The reduced efficiency of the tur- 
bulent damping of high frequency oscillations by the slow 
convective motions is taken approximately into account by 
the factor in square brackets, where we adopt s — 2. See 
|Goldrcich and Keeley (1977|), poldman and Mazeh (199l| ) 
and Goodman and Oh (1997) for a discussion. The maxi- 
mum value of the adopted viscosity coefficient v is attained 
near the lower boundary of the convective envelope. In the 
stellar models used this maximum value ranges from a few 
times 10 13 to slightly more than 10 15 in cgs units, whereby 
the latter value is reached for low forcing frequencies. Both 
the radiative and turbulent damping induce a phase lag 
of the tidal response with respect to the forcing by the 
companion and that gives rise to the tidal torques we are 
interested in. The torque- values both in between and near 
resonances are sensitive to the adopted turbulent viscosity 
law. 



3.4. The tidal torque 

The tidal torque integral corresponding with the rt-th term 
of the expansion (11) can be expressed as: 



r„(ff,n.) = < / / ImK(r)]r 4 dr 
Jo 



(19) 



where 7r results after the integration over tp, f is the con- 
stant in the tidal potential defined by Eq. (||), Im stands 
for imaginary part, and 



( 2 



J^Pf( M )X„(M)dM 



r^xi^dp 



(20) 



£ n is an overlap integral of the eigenfunctions X n (p) with 
Pf (/i) , the associated Legendre polynomial describing the 
variation of the dominant I = m = 2 tide. The tidal 
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torque is proportional to which for n = 1 and n s — > 
attains the value 48/5. 

The tidal torque follows by multiplying 7^ by the fac- 
tor m = 2. Once we have, for a given forcing frequency 
a and stellar angular rotation speed fi s , first solved the 



A„(x) eigenvalue problem and subsequently the radial os- 
cillation Eqs. (|13|-|16[), the resulting p' n (r) can be substi- 
tuted in the above integral to evaluate the tidal torque. 

4. Solution method 

4.1. Determination of angular eigenvalue X n (x) 

For a given stellar rotation frequency f2 s and forcing fre- 
quency cf the eigenvalues A„ can be determined by numer- 
ically integrating the two ordinary differential Eqs. (|J) 
and We have used a shooting method with fourth 

order Runge-Kut ta integration with variable stepsize (e.g. 
Press et al. 1992 ). To this end we rewrite the equations as 



viscous diffusion term requires three-level finite differenc- 
ing. It appears that, because of the short wavelength of 
the tidal response for low forcing frequencies a one needs 
a large number of meshpoints for a solar- type star: we have 
used 5000 radial points. The 2D oscillation code (Savonije 
and Papaloizou 1997] ) which takes the Coriolis force fully 



d/i 

dQn 



rax [i 



1 



1- M 2 



Qn, 



1 



~ A n I X n 



mx/i 



1 



Qn 



(21) 



(22) 



where the factor x = 20 s /<r expresses the importance of 
stellar rotation. 

To enable integration away from fj, = 1 it is convenient 



(1 — fi ) 2 Y n (fj,) and expand Y n (p) in a 
: c + ci (p - 1) + c 2 (u - l) 2 + . . . which 



to write X n (/j,) ■ 
power series Y n 

can be substituted in Eqs. ( J2lh and ( |22| ) to determine the 
coefficients Cfc. We use cq as a free scaling parameter and 
express coefficients ci to C4 in terms of Co, x and the un- 
known eigenvalue A n . Q n can then be expressed in terms 
of Y n as 



Qn = 



1 



X 2 fl 2 



mfi(l - x) Y n - (1 - 



M 2 ) 



dY n 
d/i 



For the even (in /1) ra — 2 solutions the boundary condi- 
tions are X n = at /1 = 1 and Q„ = at /i = 0. We can 
now integrate equations (^l|)-(p2|) from fj, = 1 — e (with 
e = 10~ 4 ) to /i = with an estimated value for the eigen- 
value A„. We iterate by adjusting A„ until the integrated 
value for Q n is sufficiently close to zero for fi = 0. 

4.2. Numerical solution of radial oscillation equations 

Once A„ (x) is determined by means of the above numerical 
method the radial oscillation Eqs. (13 16j) can be solved 
for the same value of a after prescribing the usual bound- 
ary conditions £ r = F' = at the stellar centre and by 
requiring that the Lagrangian perturbations at the stel- 
lar surface obey: 6P' = and SF'/F = 4 ST'/T (Stefan's 
law) . By transforming the differential equations into alge- 
braic equations by means of finite differences on a spatial 
mesh and applying matrix inversion similar to standard 
Henyey schemes for stellar evolution the radial equations 
can be solved (e.g. Savonije and Papaloizou 1983). The 



into account would require unrealistically long integration 
times for such a fine radial grid in combination with an 
adequate resolution in $. In order to obtain consistent 
results for the resonances we have used the traditional ap- 
proximation also outside the low frequency inertial range 
(although in Figure 2 we only show results for |a;| < 1). 

5. Stellar models 



We constructed the unperturbed ste llar models for the 
solar- type star with a recent version ( Pols et al. 1995|) o f 
the stellar evolution code developed by Egglcton (1972). 
The models represent a spherical main sequence star of 
1 Mq with chemical composition given by various val- 
ues of the central hydrogen abundance X c and Z = 0.02. 
The model was constructed by using the OPAL opacities 
( [glesias and Rogers 1996). Fig. [I] shows a frequency re- 
lated to the Brunt- Vaisala, frequency TV, namely vbv = 
sign AT I Af I in units of the stellar break-up speed O c 
versus mesh point number for both an unevolved stellar 
model and for a model near the end of core hydrogen burn- 
ing. In the latter model the square of the Brunt- Vaisala 
frequency has become quite large in the contracted helium 
core. The radiative core, where the tidal oscillations oc- 
cur, is made to contain roughly 4000 meshpoints in order 
to resolve the short wavelength response for low forcing 
frequencies. 



40 
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u 
= 
o 

cr 



-20 



-40 




vbw (*i) 

^Lamb (tl) 
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1 2 3 4 5 
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Fig. 1. The modified Brunt- Vaisala frequency ^bv = 
signA/yj^ATf and Lamb- frequency in units of the stellar 
break-up speed O c versus meshpoint number for (tx : 
0.15 Gyr) an unevolved 1 M Q model, and (t 2 : 11 Gyr) for 
a stellar model with central hydrogen abundance X c ~ 0. 
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6. Results 



6.1.3. i-mode solutions 



6.1. The solution for the eigenvalues X n 

We only show the solution of the eigenvalue problem in the 
the low-frequency inertial regime, for which the absolute 
value of the rotation parameter x = 2SI s /<j is larger than 
unity. Negative x- values correspond to retrograde, positive 
values to prograde forcing. 



6.1.1. g-mode solutions 

The calculated positive eigenvalues A„, with n — 1 to 
5, are shown as a function of 1/x in Fig. |2|a. For large 
values of \l/x\ the eigenvalues A„ approach the values 
1(1 + 1) (with l=2n), corresponding to vanishing rotation. 
The values for Ai are consistent with those calculated in 



Papaloizou and Savonije (1997). It can be seen that ro- 
tational effects become strong when \l/x\ — > 0, i.e. when 
the oscillation frequency becomes small compared to the 
angular speed of the star. 

The large values for A„ for | | 1 — > are associ- 
ated with eigenfunctions X n which attain small values for 
1 > fj, >> |x| _1 , while the peak(s) shift(s) to fj, = 0. 
This corresponds to the confinement of g-modes towards 
the stellar equator due to strong Coriolis effects near the 
poles for small oscillation frequencies. For the lowest eigen- 
value this effect is strongest for the retrograde g-modes. 
Fig. ||c shows the angular overlap integral ( |20| ) 10 log[(^] as 
a function of 1/x for A„ > 0. Due to the equatorial confine- 
ment of the g-modes their coupling to the tides becomes 
weak when \l/x\ becomes small, with strongly reduced 
values of ( n - The latter effect is enhanced due to the fact 
that the retrograde solutions have one node more in the 
•Q direction than the corresponding prograde solutions. 



6.1.2. r and q-mode solutions 

For —1/6 < x^ 1 < there are solutions with small 
eigenvalues, see Fig. ||a. The solution with A„= is a 
purely toroidal oscillation with frequency corresponding 
to x = —6. This solution corresponds to the class of r- 



modes (Papaloizou and Pringle 1978) with frequencies 
2mO B 



1(1 + 1) 

For the lowest order I = 3 (which couples to the dominant 
I = m = 2 tide) this frequency coincides with 1/x — —1/6. 
In previous papers flWS99a| ; |WS99b| ; |WS0j] ) we loosely 
called the quasi-toroidal modes (with associated small but 
non-zero values of A„) also r-modes, but here (see section 
6.5) we have identified these modes up to high radial order 
for which the non-toroidal character becomes significant 
and we denote them as q-modes. Like g-modes, the peak(s) 
of the eigenfunctions for q-modes shift to small values of 
/i when x^ 1 —* 0, although the equatorial confinement is 
less pronounced. It can be seen in Fig. ||c that the q-modes 
couple effectively to the I = m = 2 tide, with the angular 
overlap integral ( 2 ~ 10. 



Solutions with A„ < are associated with rotationally 
governed inertial i-modes which can propagate in the ap- 
proximately adiabatic convective envelope. Fig. ||b shows 
the absolute value of the eigenvalues associated with iner- 
tial solutions for m = 2 and with the adopted even symme- 
try in n for the eigenfunctions X n . When — * 1 the 
eigenfunctions of the i-solutions become more and more 
confined to the stellar poles, i.e. the coupling with the 
tide decreases when the effect of rotation becomes weaker 
(see Fig. ||d) . The eigenvalues for the prograde i-solutions 
vanish when x^ 1 — » 0, while the eigenvalues for retrograde 
i-solutions tend to zero for small, but finite values of |a;| _1 . 
The eigenvalue corresponding to the fundamental inertial 
solution drops in fact to zero for 1/x — > —1/6. At this 
point the fundamental inertial solution connects to the 
toroidal r-solution. The prograde solutions have one node 
more in the i? direction than the corresponding retrograde 
solutions. 



6.2. The eigenfrequencies: effects of stellar evolution 

By substituting the calculated values for the angular- 
eigenvalue Ai (for the given value of x = 2Q s /a) in 
Eqs. (13)-(|lrj) and solving these equations numerically 
for a range of forcing frequencies a, the resonances with 
stellar oscillation modes are readily found. By doing this 
for several stellar input models, at various phases of core 
hydrogen burning, one can follow the evolution of the nor- 
mal mode resonances of the 1 M© star from the ZAMS to 
the end of core hydrogen burning. We have calculated the 
resonances for a non-rotating star, as well as for a star 
rotating at fl s — 0.1 and 0.2 (in units of Q c ). Figs. Ha-||d 
show the evolution of the resonance frequencies for pro- 
grade and retrograde g-modes for stellar rotation speeds 
il s = 0.1 and 0.2, as the 1 M star evolves away from 
the ZAMS. In both cases the g-mode frequencies increase 
with age, as is expected from the fact that the stellar core 
contracts, so that the acceleration of gravity increases in 
the propagation region of the modes. We also find that 
the retrograde g-modes have higher resonance frequencies 
(in absolute value) than their prograde counterparts. 

Figs. ||e and ||f show the more complex evolution of the 
quasi-toroidal q-modes when the star evolves away from 
the ZAMS. The frequency of toroidal modes is propor- 
tional to the the rotation speed fl s . When the star evolves 
away from the ZAMS its radius increases, at first slowly, 
but near the end of core hydrogen burning more quickly. 
Rotation at fixed f2 s /f2 c thus means that the stellar rota- 
tion rate assumed for this figure slows in fact down with 
age because our frequency unit fl c decreases when the star 
expands. This tends to bring down the mode frequencies 
with age, as is indeed the case for the low radial order 
q-modes. However, the higher radial order q-modes cor- 
respond with larger eigenvalues Ai and thus with an in- 
creasing non-toroidal character as the radial node num- 
ber k rises, for which buoyancy effects become important. 
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Fig. 2. a The calculated positive eigenvalues A n (for n = 1 to 5) versus l/x — a/(2il s ) for m=2, corresponding to 
eigenfunctions with even symmetry in (i. The full curves correspond to n = 1, the uppermost curves to n = 5. b The 
same as panel a but now for the negative eigenvalues, c The square of the tidal angular overlap integral £ n for A„ > 
versus l/x with same line coding for n=l to 5. d The square of for A„ < versus l/x with same line coding for 
n=l to 5. 



This explains why the frequency of the high order q-mode 
resonances increase with age, like g-modes. Comparing 
Figs. |e and |f shows that the resonance frequencies of 
low-order q- modes do indeed scale with f2 s , and less so for 
the high order q-modes. 



6.3. Determination of the tidal torque 



Although the expansion flll| ) should be made over all n val- 
ues, the strength of the tidal coupling decreases (Figs. ||c- 
d) rapidly for higher orders, so that we will from now on 
consider only the first eigenvalue Ai, and ignore the intrin- 
sically weaker coupling of the higher order eigenfunctions. 

For a given forcing frequency a and stellar spin rate tt s 
we first determine Ai(a;) and substitute this in the set of 
Eqs. (p^-p^). After solving these equations (section (4.2) 
and substituting the obtained imaginary part of p' in the 
torque integral ( |l^ ) the tidal torque can be evaluated. The 
calculations are done for a perturbing companion of 1 M 
in circular orbit with binary separation D = 4i? s , i.e. for 
a fixed ratio of R s /D. The values for the torque integral 



71 have to be multiplied by a factor m — 2 to get the tidal 
torque T. 

The tidal torque on a non-synchronously rotating 
solar-type star in a given binary system with a compan- 
ion of mass M p in a circular orbit with arbitrary orbital 
separation D is then found from 



T(cr) = m 



Mp 



2 / 4j?g X 6 

D 



where, for a stellar model with X c ~ 0.4 and a given value 
of (j — (j — mfl s with m = 2, T\ can be read from Figs. |J 
to 0. 



6.4. Comparison with \T98^ s calculation 

To check our results we have compared some of our torque- 



values for the non-rotating case with those of T98 



there exist no detailed results for rotating stars to com- 
pare with. Table gives some results of both calculations 
for three different circular orbits with periods of 4.2, 8.5 
and 12.4 days. In the table we have scaled our results to 
those of T98, i.e. we use their normalisation of the tidal 
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Fig. 3. The variation of resonance frequencies Co corresponding to Ai > with stellar evolution, a and b The prograde 
g-modes for rotation rate f2 s = 0.1 and 0.2, respectively, c and d The retrograde g-modes for rotation rate O s = 0.1 
and 0.2. e and f The quasi-toroidal q-modes for rotation rate Q s = 0.1 and 0.2. 



torques as in their Fig. 2. After comparing their torque 
values with our results for non-rotating stars we conclude 
that the two calculations yield quite similar results. 

6.5. Torque spectra for Ai > (g- and q-modes) 

After repeating the torque calculation procedure for a 
wide range of frequencies (\a\ < 1, where all frequencies 
stated without units arc in units of the critical rotation 



rate f2 c ), we end up with a number of torque spectra, like 
the one shown in Fig. [| for fl s = 0.1 and X c = 0.4. Note 
that the actual height of the resonances depends on the 
adopted turbulent viscosity law for the convective regions 
and on the amount of turbulent overshooting. We have 
assumed a Gaussian decay of the viscosity coefficient be- 
yond the convection boundaries with a decay length of 0.5 
pressure scale height. 
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Table 1. Comparison of calculated torque values: the up- 
per values are from T9S , while the lower values correspond 
to the calculations presented here. T pca k and T a g are the 
torque values at the resonance peaks considered, and for 
neighbouring off-resonant frequencies, respectively. The 
middle column gives the relative FWHM of the resonance 
peaks. 



P(d) 


Tpoak (Cgs) 


Aa/a 


Toft (cgs) 


4.2 


3.2 x 10 38 
3.4 x 10 38 


3 x 10" 6 
7 x 10" 6 


4 x 10 35 
7 x 10 35 


8.5 


2.8 x 10 35 
3.6 x 10 35 


4 x 10" 5 
4 x 10" 5 


6 x 10 34 
5 x 10 34 


12.4 


6.3 x 10 33 
1.5 x 10 34 


2 x 10" 4 
2 x 10" 4 


5 x 10 33 
7 x 10 33 



The high frequencies near \a\ — 1 will never occur 
in normal binaries with solar-type stars, but part of the 
higher frequencies may be relevant to extreme systems 
with planets in very eccentric orbits. For \a\ < 0.25 the 
resonance spectrum of prograde and retrograde high radial 
order g-modes becomes so dense that in Fig. || it can no 
longer be resolved. For decreasing forcing frequency \a\ 
the radial wavelength of the g-modes decreases, so that 
the resonances fade away by severe radiative damping. 
Also the tidal coupling strength £i of (especially the ret- 
rograde) g-modes declines with \a\ (see Fig. For the 
high frequency g-mode resonances turbulent viscosity in 
the (lower part) of the convective envelope provides the 
dominant damping. The turnover to predominant radia- 
tive damping occurs near \a\ ~ 0.1 

We have zoomed in to determine the resonances with 
the stellar normal modes down to the (somewhat arbi- 
trary) frequency where the resonances become too weak 
to trace down individually. 

For the prograde g-modes the resonance peaks beyond 
about g4so at u — 0.03 (1/x ~ 0.15) have been ignored, 
while on the retrograde side of corotation g-modes with 
up to about 800 radial nodes have been identified. 

The off-resonant torque values follow a curved line 
which reflects two main effects: for high oscillation fre- 
quencies \a\ the turbulent dissipation in the convective 
envelope becomes less efficient due to the mismatch of the 
low frequency tur bule nt motions and the high frequency 
oscillations (Sect. |3.3| ), reducing the off-resonant torque- 
values. On the other hand, for low oscillation frequencies 
\a\ < 0.1 the peak torque- values decrease strongly by the 
combined effect of enhanced radiative damping and the 
rotational confinement of the g-mode response to regions 
close to the stellar equator. 

When we increase the (uncertain) turbulent viscosity 
in the convective envelope by an order of magnitude the 
off-resonant torque values increase by roughly an order of 
magnitude, while the peak values at the resonances (for 
\a\ > 0.1 where radiative damping is weaker) decrease by a 



similar factor. For low forcing frequencies the tidal g-mode 
response has short radial wavelength, and the resonances 
are predominantly damped by radiative damping (for the 
adopted turbulent viscosity law). 

For fi s = 0.1 the torque distribution (Fig. [|) shows 
a dip near a ~ —0.11. The same dip occurs more promi- 
nently at about the same value of 1/x ~ —0.5 for il s = 0.2 
(or higher rates) , see Fig. When the convective envelope 
is artificially removed from the star the dip in the torque 
distribution disappears. This indicates the dip is related 
to rotation and is caused by the convective envelope. 

For frequencies — f2 s /3 < a < the eigenvalues Ai are 
small (see Fig. ||a) . These small eigenvalues are associated 
with the closely spaced quasi-toroidal q-modes which can 
be ordered as q^, with k denoting the number of radial 
nodes. Because of the long wavelength of the lower ra- 
dial order q-modes these retrograde oscillations are weakly 
damped. They couple efficiently to the I = m = 2 dom- 
inant tidal component: the angular overlap integral £i is 
relatively large, except when — > 0. In Fig. [I] and || 

the ~ 530 identified q-mode resonances are not individ- 
ually resolved and produce the black peak situated just 
left of corotation (a = 0). Although the damping of es- 
pecially the low radial order (with small A) oscillations of 
this type may be underestimated by our simple treatment, 
the several hundred relatively weakly damped resonances 
with q-modes provide a significant tidal effect by which 
a rapidly spinning solar-type star can be spun down effi- 
ciently. 

6.6. Torque spectra for Ai < (inertia I modes) 

For negative values of the eigenvalue Ai the corresponding 
tidal oscillations are restricted to the convective envelope 
(in contrast to all other oscillation modes considered here, 
which are all limited to the radiative stellar core) and 
are associated with a spectrum of rotationally governed 
inertial modes. These oscillation modes with frequencies 
limited to the inertial range (\a\ < 2f2 s ) produce an ex- 
tended, relatively high torque level including a strong res- 
onance peak (broadened by the turbulent viscosity) in the 
frequency range where the torque due to retrograde g- 
modes falls off sharply, filling the gap left of the q-modes: 
compare Figs. ^ and [|. The height of the inertial torque 
peak increases with the stellar rotation rate O s and de- 
creases when the turbulent viscosity is increased. Outside 
the strong resonance peak the torque- values increase with 
increasing viscosity. For the low-frequencies in the inertial 
range the value of the reduction factor in the expression 
for the turbulent viscosity coefficient Eq. ([l8]) is practically 
equal to unity. The relatively strong inertial response can 
be an important driving factor for the establishment of 
resonance locking (|WS99a|; |WS99b|; |WS0l|). 



7. Conclusions 

We have calculated the tidal torque on a rotating 1 Mq 
star due to an orbiting companion, using the traditional 
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Fig. 4. The torque integral T versus forcing frequency a for Ai > 0. The calculation is for a 1 Mq star with central 
hydrogen abundance X c — 0.4, rotating at speed fi s = 0.1 f2 c - The narrow dip before each (low order) resonance peak 
has been ignored in this plot. Negative values of a correspond to retrograde forcing, for which the sign of the torque 
should be read as negative. Results are for a fixed ratio of orbital separation to radius D / R s = 4 and companion mass 



M„ 



1 Mq. Many strong resonances with quasi-toroidal q- modes, shown here as an unresolved black peak just left 



of (7 = 0, can be identified (with up to ~ 500 radial nodes) for this stellar model. 
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Fig. 5. The torque integral versus forcing frequency for Ai > 0. Same as Fig. Q but now for rotation rate fl s = 0.2£l c 
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Fig. 6. The tidal torque integral versus forcing frequency 
for Ai < 0, (i.e. due to excitation of inertial modes in 
the convective envelope) for several stellar rotation speeds: 
n s = 0.05, 0.10, 0.15 and 0.20 in units of fi c . The dashed 
lines correspond to f2 s =0.05, the dot-dashed lines to 0.10 
etc. The companion mass is 1 M Q , while the ratio D/R s = 
4 was again fixed. 

approximation. We have been able to identify both g-mode 
and quasi-toroidal q-mode resonances with up to ~ 1000 
radial nodes in the more evolved main sequence models by 
using 5000 meshpoints for the radial grid and calculated 
the corresponding tidal torque spectra for stellar models 
up to the end of core hydrogen burning. The tidal torques 
calculated in this paper are based on a standard local 
mixing length approximation for turbulent convection in 
the convective envelope of solar-type stars, with a simple 
prescription for reduced viscous damping at high forcing 
frequencies. The resonances with low radial order tidal 
oscillations are damped predominantly by the turbulent 
viscosity in the lower convective envelope, and the cor- 
responding peak values for the tidal torque are inversely 
proportional to the adopted coefficient of turbulent vis- 
cosity. Radiative diffusion in the radial direction is also 
contributing to the tidal torque (strong radiative diffusion 
dominates the damping of resonances with high radial or- 
der g-modes at low forcing frequencies which results in 
weak tidal torques). 

For low retrograde forcing frequencies we find a rel- 
atively strong tidal response in the convective envelope 
for angular eigenvalues Ai < 0, which corresponds to tur- 
bulent dissipation of tidally excited inertial modes. This 
spectrum of inertial modes may provide an efficient driv- 
ing mechanism for tidal resonance locking and may gen- 
erate an enhanced tidal evolution in low-mass binary sys- 
tems. 

For still smaller \a\ resonances with the retrograde 
quasi-toroidal q-modes (propagating in the radiative core) 
become a prominent feature in the torque spectrum. 
This q-mode spectrum consists of several hundred strong, 
closely spaced resonances. When the binary parameters 
and stellar rotation rate are such that a low orbital har- 



monic falls in the frequency regime of the q-modes, en- 
hanced tidal evolution must occur when the harmonic 
evolves through the many closely spaced q-mode reso- 
nances. 

The effects mentioned above are related to rotation of 
the binary components and have been ignored in previous 
studies of tidal effects in low-mass binary systems. In a 
forthcoming paper we will apply the above results to study 
the tidal evolution of low-mass eccentric binary systems. 

References 

Claret, A. and Cunha, N. C. S.: 1997, A&A 318, 187 
Cowling, T. G.: 1941, MNRAS 101, 367 
Duquennoy, A. and Mayor, M.: 1991, A&A 248, 485 
Eggleton, P. P.: 1972, MNRAS 156, 361 
Goldman, I. and Mazeh, T.: 1991, ApJ 376, 260 
Goldreich, P. and Keeley, D. A.: 1977, ApJ 211, 934 
Goodman, J. and Dickson, E. S.: 1998, ApJ 507, 938 
Goodman, J. and Oh, S. P.: 1997, ApJ 486, 403 
Iglesias, C. A. and Rogers, F. J.: 1996, ApJ 464, 943 
Latham, D. W., Mazeh, T., Stefanik, R. P., Davis, R. J., 
Carney, B. W., Krymolowski, Y., Laird, J. B., Torres, 

G. , and Morse, J. A.: 1992, AJ 104, 774 

Mathicu, R. D., Duquennoy, A., Latham, D. W., Mayor, 
M., Mermilliod, T., and Mazeh, J. C: 1992, in Binaries 
as tracers of stellar formation. A. Duquennoy and M. 
Mayor (eds), p. 278, Cambridge University Press 
Papaloizou, J. and Pringlc, J. E.: 1978, MNRAS 182, 423 
Papaloizou, J. C. B. and Savonije, G. J.: 1997, MNRAS 
291, 651 

Pols, O. R., Tout, C. A., Eggleton, P. P., and Han, Z.: 

1995, MNRAS 274, 964 
Press, W. H., Tcukolsky, S. A., Vetterling, W. T., and 

Flannery, B. P.: 1992, Numerical recipes in FORTRAN. 

The art of scientific computing, 2nd ed., Cambridge 

University Press 
Savonije, G. J. and Papaloizou, J. C. B.: 1983, MNRAS 

203, 581 

Savonije, G. J. and Papaloizou, J. C. B.: 1997, MNRAS 
291, 633 

Savonije, G. J., Papaloizou, J. C. B., and Alberts, F.: 1995, 

MNRAS 277, 471 
Terquem, C, Papaloizou, J. C. B., Nelson, R. P., and Lin, 

D. N. C: 1998, ApJ 502, 788 
Unno, W., Osaki, Y., Ando, H., Saio, H., and Shibahashi, 

H. : 1989, Nonradial oscillations of stars, Univ. of Tokyo 
Press 

Verbunt, F. and Phinncy, E. S.: 1995, A&A 296, 709 
Witte, M. G. and Savonije, G. J.: 1999a, A&A 341, 842 
Witte, M. G. and Savonije, G. J.: 1999b, A&A 350, 129 
Witte, M. G. and Savonije, G. J.: 2001, A&A 366, 840 
Zahn, J. P.: 1977, A&A 57, 383 



